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lO ■ Abstract 



The purpose of this paper is to study quasinormal modes (QNM) of the Bardeen 
black hole due to scalar perturbations. We have done a thorough analysis of the 
QNM frequencies by varying the charge q, mass M and the spherical harmonic index 
/. The unstable null geodesies are used to compute the QNM's in the eikonal limit. 
Furthermore, massive scalar field modes are also studied by varying the mass of the 
{Sj \ field. Comparisons are done with the QNM frequencies of the Reissner- Nordstrom 

black hole. 
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1 Introduction 

In general, black hole space-times are expected to have horizons as well as singularities 
^ ■ covered by the horizons. Contrary to this notion, a "regular" space-time without a 

singularity and with a horizon was proposed by Bardeen [1]. This particular paper 
is not readily available. However, a discussion of this model is given by Borde in [2] 
[3]. There were other regular black holes other than the one proposed by Bardeen 
and they were all referred to as " Bardeen black holes" by Borde [2]. In this paper, 
we will focus on the space-time proposed by Bardeen [1]. 

Ayon-Beato and Garcia [4] proposed a model of nonlinear electrodynamics cou- 
pled to Einstein gravity to obtain Bardeen black hole as an exact solution. Hence, 
the Bardeen black hole can be interpreted as the solution to a nonlinear magnetic 
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monopole with a mass M and a charge q [4]. Ayon-Beato and Garcia [4] have pre- 
sented several other interesting regular solutions with nonlinear electrodynamics cou- 
pled to General Relativity in [5] [6] [7] [8]. 

There are several works in the literature related to the Bardeen black hole. Grav- 
itational lensing of the regular black hole was studied by Eiroa and Sendra [9]. The 
geodesic structure of the test particles around the Bardeen black hole were studied by 
Zhou et.al [10]. Gravitational and electromagnetic stability were discussed by Moreno 
and Sarbach [11]. Quantum corrections for the Bardeen black hole was presented by 
Sharif and Javed [12]. 

In this paper, our focus is on studying the scalar field perturbations of the Bardeen 
black hole and to compute the quasinormal modes of the perturbations. 

When a black hole undergo perturbations, the resulting behavior can be described 
in three stages. The first stage corresponds to radiation due to the initial conditions of 
the perturbations. The second stage corresponds to damped oscillations with complex 
frequencies. These frequencies are independent of the initial conditions and are only 
dependent on the black hole properties such as the mass, charge and the angular 
momentum. These modes are called quasinormal modes (QNM). The third stage in 
general corresponds to a power law decay of the fields. 

QNM's have attracted lot of attention from the research community. First of 
all, there is interest from the experimental point of view, since, there is hope that 
the QNM's may be detected by the gravitational antennas such as LIGO, VIRGO 
and LISA in the future. Since the QNM's only depend on the properties of the 
black holes, such detections would give clues to identify the physical properties of 
the black holes. On the other hand, due to the famous relation between AdS/CFT 
duality, many works have focused on studying QNM's of black holes with a negative 
cosmological constant [13] [14]. Another reason to create interest on QNM's was the 
conjecture by Hod [15] relating quantum properties of Schwarzschild black hole and 
asymptotic QNM's. There are many works aimed at computing asymptotic QNM 
frequencies along those lines. An excellent review on QNM's is written by Konoplya 
and Zhidenko [16]. It is fascinating to see some new work relating QNM to various 
aspects of black hole physics. For example, some recent work have addressed relation 
between QNM's and hidden conformal symmetry [17] [18]. Another interesting paper 
was written on the connection between gravitational lensing and QNM's by Stefanov 
et.al [19]. Not only black holes, even the naked singularities have been studied from 
the QNM point of view [20]. Given all the above, it is worthwhile to seek answers 
how nonlinear sources modify the QNM properties of a black hole. 

The paper is presented as follows: In section 2, the Bardeen black hole solutions 
are introduced. In section 3, the perturbations by a massless scalar filed is given. 
In section 4, we will computer the QNM's using the sixth order WKB approach and 
discuss the results. In section 5, a relation between the null geodesies and the QNM's 
are presented. In section 6, the perturbations by a massive scalar field is studied. 
The summary is given in section 7. Directions for further studies are given in section 
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2 Introduction to the regular Bardeen black hole 

In this section, we will give an introduction to the regular static charged black hole 
named as Bardeen black hole [1]. Ayon-Beato and Garcia [4] interpreted the Bardeen 
black hole as the gravitational field of a magnetic monopole arising from non-linear 
electrodynamics. The proposed action to include the non-linear electrodynamic term 
is, 

s =/ d4 ^[i^G-s £ H (1 > 

Here, R is the scalar curvature, and C(F) is a function of F — ^F^ U F^ U . Here, 
F^v = 2(\7^A V — XJpA^) is the electromagnetic field strength. In [4] the authors 
derived the function C(F) in order to obtain the Bardeen black hole as, 



Here, q and M are the magnetic charge and the mass of the magnetic monopole. 
Also, s = 4jL. The equations of motion derived from the action in eq.(l) is given by, 

g; = 2 {^F,,^ - s;cj (3) 

V, = (*) 

Static spherically symmetric solution for the above equations were proved to be the 
Bardeen black hole solution given by the metric, 

ds 2 = -f(r)dt 2 + f(r)- l dr 2 + r 2 (d6 2 + sin 2 (6)dip 2 ) (5) 

where, 

2Mr 2 

/(r) = 1 " p + g2)3 /2 (6) 

The magnetic field is given by, 

F dv = 2qsin9 (7) 

For q 7^ 0, the space-time in eq.(5) has horizons only if, |g| < |^|. This was shown by 
Borde [2] [3]. For q > |^=, there are no horizons. For q = |^=, there are degenerate 
horizons. The function f(r) is plotted in Fig.l for varying magnetic charge q. 
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Figure 1. The figure shows the function f(r) for M = 1 and varying values of q. 

Asymptotically, the metric function /(r) behaves as, 

2M 3Mq 2 ^ / 1 \ 

The metric for the black hole in Einstein-Maxwell gravity, given by the well known 
Reissner-Nordstrom black hole, with a magnetic charge is, 

ds 2 = -f{r) RN dt 2 + f(r)~ R l N dr 2 + r 2 {d6 2 + sin 2 (6)d V 2 ) (9) 

where, 

/(r)RJV = 1 _^ + £! (10) 

In Fig.2, the two metric functions for the Bardeen black hole and the Reissner- 
Nordstrom black hole are plotted for comparison. It is clear from the Fig.2, that 
both black holes have two horizons. For small r, the behavior is some what different 
even though asymptotically, both functions are similar. The non-singular nature of 
the function f(r) for the Bardeen black hole is observed from the Fig.2. As dis- 
cussed in [4], the space-time is regular everywhere since all the scalar curvatures, 
R, R I1U R^ U and R l t VOl pR liva ^ are regular every where. However, the electromagnetic 

2 

invariant F = ^ has singular behavior. The Bardeen black hole satisfy the weak 
energy condition. 
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Figure 2. The figure shows the function f(r) for the Bardeen black hole(dark) 
and the Reissner- Nordstrom black hole(dashed). Here, M = 1 and q = 0.5. 



The Hawking temperature of the Bardeen black hole is given by, 

'2Mr + (r 2 , - 2q 2 ) 



rp l _ d 9tt 



An dr 



1 

4tt 



(q 2 + 4)5/2) 



11' 



Here, r + is the event horizon of the black hole which is a solution of f(r) = 0. In 
comparison, the Hawking temperature for the Reissner-Nordstrom black hole is, 



Trn ~ 



2M 



r 3 



(12) 



The temperature for both black holes for the same mass is plotted in Fig.3. The 
Reissner-Nordstrom black hole is "hotter" than the Bardeen black hole. 
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Figure 3. The figure shows the temperature T for the Bardeen black hole (dark) 
and the Reissner-Nordstrom black hole (dashed) as a function of the magnetic charge 
q. Here, M = 1. 



3 Masslesss scalar perturbation of Bardeen black 
holes 

In this section, we will introduce scalar perturbation by a massless field around the 
Bardeen black hole. The Klein-Gordon equation for a massless scalar field $ in curved 
space-time can be written as, 

V 2 $ = (13) 

which is equal to, 

-^d^V^gg^®) = (14) 



Using the ansatz for the scalar field <3>, 

$ = e -^y(0,0)^ (15) 
r 

eq.(14) simplifies to the Schrodinger-type equation given by, 

drl 

where, 



+ (u 2 - V(n)) £(r) = (16) 



V(r) = Kl + y (r) + /(r)f(r) (17) 
Here, r* is the well known "tortoise" coordinate given by, 

dr , . 

* = IF) (18) 

Note that / is the spherical harmonic index. Here, r* cannot be evaluated explicitly 
due to the nature fo the function f(r). When r — > oo, r* — > oo and when r — > r + , 
r* — > —oo. 

The effective potential V(r) for the Bardeen black hole is plotted to display how 
it changes with charge q, the mass M, and the spherical harmonic index / in Fig. 4, 
Fig. 5 and Fig. 6 respectively. 
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Figure 4. The behavior of the effective potential V(r) with the charge for the 
Bardeen black hole. Here, M — 1 and 1 = 2. The height of the potential decreases 
when the charge decreases. 
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Figure 5. The behavior of the effective potential V(r) with the mass M. Here, 
q = 0.5 and I = 2. The maximum height of the potential increases as M decreases. 
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Figure 6. The behavior of the effective potential V(r) with the spherical harmonic 
index I. Here, M = 1 and q = 0.5. The height of the potential increases when / 
increases. 

We have also plotted the scalar effective potential for the Reissner-Nordstrom 
black hole with the one for the Bardeen black hole in Fig. 7 for comparison. 
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Figure 7. The behavior of the effective potential V(r) for the Bardeen black hole 
(dark) in comparison with the potential for the Reissner-Nordstrom black hole(dashed). 
Here, M = 1, q = 0.2 and I = 2. The potential for the Bardeen black hole is slightly 
higher in the small r range. 

3.1 Remarks on the stability of the black hole 

The potentials are real and positive outside the event horizon for all the figures. 
Hence, following the arguments by Chandrasekhar [21] the Bardeen black holes can 
be considered stable classically under perturbations by a massless scalar field. 
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4 Quasi-normal modes of Bardeen black hole 



Quasi-normal modes (QNM) for a perturbed black hole space-times are the solutions 
to the wave equation given in eq.(16). In order to obtain solutions, one has to impose 
boundary conditions. At the horizons, the boundary condition is such that the wave 
has to be purely ingoing. In asymptotically flat space-times, such as the Bardeen 
space-time, the second boundary condition is for the solution to be purely outgoing 
at spatial infinity. 

Usually, the wave equation for black hole perturbations cannot be solved exactly. 
There are few cases of exactly solved models known to the authors which are men- 
tioned here. In 2+1 dimensions, the wave equations of the well known BTZ black 
hole [22], the charged dilaton black hole [23] [24], the Lifshitz black hole [25] and the 
Godel black hole [26] has been solved to obtain exact QNM values. In two dimen- 
sions, an asymptotically anti-de-Sitter black hole has been solved exactly [27]. In five 
dimensions, Nunez and Starinets have obtained exact values for vector perturbations 
[28]. 

There are many methods developed to compute QNM's in the literature. In this 
paper, a semi analytical technique developed by Iyer and Will [29] is followed. The 
method is based on the WKB approximation. Iyer and Will developed it up to third 
order and later, Konoplya developed it up to sixth order [30]. Konoplya computed 
QNM frequencies of D dimensional Schwarzschild black holes in that paper which also 
includes a comparison WKB method with varying orders. Examples of the application 
of the third order WKB to compute QNM' are given in [31] [32] and of the sixth order 
WKB is given in [33]. 

We will follow the formalism presented in the paper by Konoplya [30]. In the 
WKB formula, the QNM frequencies are related to the effective potential in eq.(17) 
as, 

° - L 2 - L 3 - L A - L 5 - L 6 = n + - (19) 

2 

Here, V and V " are the maximum potential and the second derivative of the potential 
where the maximum occurs. The expressions for L 2 and L 3 are given in [29] and L 3 , 
L 4 , L 5 and L 6 are given in [30]. Here, n is the overtone number. The computed u 
values are complex and are as oj = ojr — iuuj. 

First, we have computed QNM frequencies by varying the charge of the black hole. 
We have also computed the QNM's for the Reissner-Nordstrom black hole with the 
same mass and the charge in order to compare. Here, n — 0. 
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Figure 8. The behavior of Re u with the magnetic charge q for M = 1, and / 
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Figure 9. The behavior of Im u with the magnetic charge q for M = 1, and I 
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From the above Fig. 8 and Fig. 9, it is clear that real value of the QNM frequency 
oj increases when q increases for both black holes. However, the imaginary part of u 
decreases for Bardeen black hole with charge, while for the Reissner-Nordstrom black 
hole, there is a maximum before it starts to decrease. We like to mention here that 
QNM frequencies for Reissner-Nordstrom black hole were studied by Anderson [34] 
and Leaver [35]. There is a discussion on some of the earlier work on this subject in 
the book by Frolov and Novikov [36]. 

Next, the QNM values are computed for the Bardeen black hole for various values 
of the mass M as shown in the following figures, Fig. 10 and Fig. 11. Here, n = 0. 
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Figure 10. The behavior of Re u with the mass M for I = 2, and q = 0.5. 
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Figure 11. The behavior of Im oj with the mass M for 1 = 2, and q = 0.5. 
When mass M is increases, both ur and ui decreases. 

Next, the QNM values are computed for the Bardeen black hole for various values 
of the spherical harmonic index / as shown in the following figures, Fig. 12, Fig. 13, 
Fig. 14 and Fig. 15. Note that the WKB work only for I > n. Hence, we have chosen 
I = 2 and computed u of the fundamental mode with n = and the first over tone 
with n = 1. 
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Figure 12. The behavior of Re oo with the spherical harmonic index I for M — 1, 
and q = 0.5. Here n = 



0.0940 




Figure 13. The behavior of Im oo with the spherical harmonic index I for M = 1, 
and g = 0.5. Here n = 0. 
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Figure 14. The behavior of Re u with the spherical harmonic index I for M — 1, 
and q = 0.5. here n = 1. 
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Figure 15. The behavior of Im u with the spherical harmonic index / for M = 1, 
and q = 0.5. here n = 1. 

For both n = and n = 1, ujr increases linearly with /. On the other hand, uj 
decreases and becomes stable for large I. 

5 QNM's of the massless black holes in the eikonal 
limit from the null geodesies of the black hole 

5.1 uj to the lowest order via null geodesies 

Cardoso et.al [37] presented an important result to compute the QNM frequencies at 
the eikonal limit via the unstable null geodesies of the black hole for asymptotically 
flat black holes. This method was based on some earlier work done along these lines 
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by Mashhon et.al[38] [39]. This approach has been applied to the Kerr black hole 
by Dolan[40], near extreme Kerr black hole by Hod [44] and to the black holes in 
anti-de-Sitter space by Morgan et.al [45]. 

First, let us give an introduction to the null geodesies of the Bardeen black hole. 
The geodesies of the Bardeen black hole were studied in detail by Zhou et.al [10]. 
Hence referring to further details to that paper, we will only present the final equation 
of motion of the photons as, 

r 2 + V null = E 2 (20) 

with, 

Vnull = (*pj f(r) (21) 

Here, L is the angular momentum of the photons. For r = r^, V nu ii = and for 
r — > oo, V nu u — > 0. In the Fig. 13, the V nu ii is given for various values of the magnetic 
charge q. The height is higher for higher charge q. 



0.115 




Figure 16. The graph shows the relation of V nu ii with r for various values of the 
magnetic charge q. Here, M = 1 and L = 1. When q increases, Vnuii increases. 

Referring to the effective potential for the massless scalar field V SC aiar in eq.(17), 
one can see that in the eikonal limit (I — > oo), 

Vgcalar ~ ~~ ^ (22) 

Hence, one can conclude that the maximum of V sca i ar occurs at r = r m given by, 

2/(r TO ) - r m f{r m ) = (23) 

T 2 f 

Since the effective potential for the null geodesies, is given by V nu u = -pf, the maxi- 
mum of V nu u occurs at Vn ull = leading to, 

2/(r c ) - r c f(r c ) = (24) 
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Hence the maximum of V sca i ar and the location of the maximum of the null geodesies 
coincides at r m = r c . Cardoso et.al. [37] presented the QNM frequencies in the eikonal 
limit, as, 

uqnm = tt c l - i(n + -)\\\ (25) 



Here, fl c is the coordinate angular velocity given as, 



n c = 4 
t 



(26) 



and, A is the Lyapunov exponent which is interpreted as the decay rate of the unstable 
circular null geodesies. The derivation of the above results is clearly given in Cardoso 
et.al[37]. For the Bardeen black hole, Q c and A are given as, 



f(r c 



i{r c ) \ rl \ r2( r 2 + g2)3/2 



( r 2 + g2 ) 3 /2 _ 2Mr c 2 



(27) 



A 



\l 2i(r c )2 



2L 2 



(2f 



From the values of Q c and A, one can extract the real and the imaginary part of u 
from eq.(25) easily. In the Fig. 17 and Fig. 18, Q and A are given for the Bardeen 
black hole. One can conclude that the behavior is very similar when u was computed 
using the WKB approach in section 4. 
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Figure 17. The graph shows Vt c as a function of q. Here, M = 1 and L = 1 
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Figure 18. The graph shows the Lyapunov exponent A as a function of q. Here, 
M = 1 and L = 1. 

One can compare the oj in the eikonal limit for the Reissner-Nordstrom black hole 
with the Bardeen black hole. First let us compute the r c for the maximum of the 
effective potential for the photons of the Reissner-Nordstrom black hole as, 



Hence, Xrn and Qrn can be computed with the eq.(27) and eq.(28) with the corre- 
sponding metric function /W( r )- In the Fig.19 and Fig.20 \rn and f2 RA r are plotted. 
Note that \rn has a maximum at q — 0.7 M. For the Reissner-Nordstrom black hole, 
the oscillation frequency, ur, increases with the charge q. The damping rate, ui, 
increases and reaches a maximum at q — 0.7 M to decrease rapidly to zero at M — q. 
Ferrari and Mashhoon [39] presented QNM's of the Reissner-Nordstrom black hole. In 
that paper, the QNM's were presented in the eikonal limit ( I 3> 1) and showed similar 
behavior. We like to note that the Lyapunov coefficients for the Reissner-Nordstrom 
black hole has been studied in detail in a recent paper by Pradhan [46]. 




(29) 
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Figure 19. The graph shows f2 c as a function of q for the Reissner-Nordstrom black 
hole. Here, M — 1 and L = 1 
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Figure 20. The graph shows the Lyapunov exponent A as a function of q for the 
Reissner-Nordstrom black hole. Here, M = 1 and L — 1. 



5.2 On the expansion method to compute the QNM's at the 
eikonal limit 

Dolan and Ottewill [41] developed an expansion method to compute QNM's at the 
eikonal limit to higher orders. Here, the basis for their approach is basically the same 
as given in section( 5.1) where the null geodesies and the unstable circular orbits are 
the key to the computation. Other references where this method is applied are given 
in [42] [43] [40]. 

In this approach, the QNM's were expanded in inverse powers of L = I + 1/2. Let 
us first outline the formalism so that later we can apply it to the Bardeen black hole. 
We shall start with the master equation for the scalar perturbation given in eq.(16), 
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drl 

where, 



+ (co 2 - V scalar (n)) £(r) = (30) 



Vscalar{r) = — 1 (31) 

Now, redefine £(r) as, 

£(r) = expW 7$*) (32) 

Here, 

a(r) = ib c k c (r)u (33) 

The parameter b is called the "impact parameter" given by L/E and b c is the value at 
r = r c which is the radius of the unstable circular orbit. L and E of the null geodesies 
are conserved quantities of the orbits given by, 

L = r 2 0; E = f(r)t (34) 

More information about the geodesies are given in [10] [41] if one needs more details 
of the origin of these quantities. 

Now, to define what k c (r) is, a new function k 2 (r,b) is defined as, 

k\r,b) = l 2 -{^± (35) 



r 



The origin of the function is at in the equation for null geodesies given in eq.(20) 
which also could be written as, 

\d(j>) r 2 b 2 r 2 V } 

In Dolan and Ottewill's paper [41], the assumption is made that, there exists a critical 
impact paramter b c , such that k 2 (r,b c ) has degenerated roots. Since at r = r c , 
dr/d(j)= 0, this leads to, 



dr 



= (37) 



r=rc 



If the assumption is made that the repeated root is a double root, a new function 
k c (r) is defined as, 

k c (r) = Sign{r - r c )sjk 2 {r,b c ) = (r - r c )K{r) (38) 
Here, K(r) becomes, 

K(r) 



\ (r - r c f 
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Hence, the function k c (r) is positive for r > r c and negative for r < r c . If just 
J k 2 (r, b c ) is considered to be k c (r), then it would have lead to a function which is 
not differentiable at r = r c . On the other hand, the definition given in eq.(38) leads 
to a smooth function at r = r c . To clarify these issues, we have plotted the graphs 
k 2 (r,b c ), ^Jk 2 (r, b c ), k c {r) and dk c {r)/dr in Fig. 21 and Fig. 22. 
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Figure 21. The behavior of k(r,b c ) and yk(r,b c ) 2 with r. Here M — 1, q — 0.76, 
r c = 2.3299 and b c = 4.5484. 
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Figure 22. The behavior of k c (r) and with r. Here M = 1, g = 0.76, 

r c = 2.3299 and b c = 4.5484. 

Now that k c (r) is well established, on can substitute £(r) into eq.(30) which sim- 
plifies to be, 

ar \ ar / ar 



( w 2 + a ( r )3 _ K catar (r) + /(r)a'(r)) ^ = ( 40 ) 
Following the approach in [41], u and v (r) are expanded in inverse powers of L as, 

w n =o = ^o-i + «o + L^ai + L~ 2 a 2 + (41) 



and 



;„=o = exp (S {r) + L- 1 S 1 (r) + L- 2 S 2 (r) + ....) 



(42) 



Note that in this paper, we are only interested in the fundamental modes correspond- 
ing to n — 0. The expansion could be done for higher modes with n > as given 
in [41]. Now, the u> and v (r) are substituted to eq.(40) and similar powers of L are 
collected together as, 

1 

ai 



\ 



r 2 c b c 



(43) 



L 1 : 



2a_iao 2b 2 c k c {r) 2 a-idQ 



fir) f(r) 

ib c a-ik' c (r) + 2ib c k c (r)a-iS' (r) = 



(44) 
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^T2-— +{ a o + 2a-iOi j ^ 



(a 5'i(r) + a_i5i(r))(2i6 c A; c (r))+ 

iaoM&O + f(r)^(r) + /(r) (^(rf + Sj'(r)) = (45) 

Here ' denotes differentiation with respect to r. The coefficients <2j are found by 
imposing continuity condition at r = r c . Once a, are found, S'^r) could be found. 
For example, from L 2 terms, a_i is evaluated at r = r c as, 



a_i = 



/(rc) (46) 



Note that this is same as Q c computed in eq.(27). 
From L 2 term, ao is evaluated at r = r c as, 

«o = ^f(r c )K(r c ) (47) 

By substituting a and a_i back into eq.(44), one obtains S' (r). One can differentiate 
S' (r) to obtain S'o(r) and substitute to eq.(45). Then, a\ can be evaluated r = r c . 
This process can be continued to find all Oj values and all ^'(r) functions. In this 
paper, we will only compute a_i, a and a\ for the Bardeen black hole. We have 
chosen, M = 1 and q = 0.76 leading to r c = 2.3299 and b c = 4.5484. By the 
expansion method, the following a, values are computed. 

a_i = 0.219858; a = -0.0792063i; ai = 0.00499392 (48) 
Finally, o; n= o in powers of L is written as, 

co n=0 = 0.219858L - 0.0792063i + 0.00499392ZT 1 (49) 



6 Massive scalar perturbations 

In this section, we will address how the massive scalar field decay. For the Schwarzschild 
black hole, it has been observed that the massive modes decay slower than the mass- 
less field [47] [48]. Hence it is interesting to see if such behavior is possible in the 
Bardeen black hole, leading to long lived modes. 

Let us first present the equation for a massive scalar field in curved space-time as, 

V 2 $ - m 2 $ = (50) 
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Using the ansatz similar to eq.(15), one obtain the modified potential as, 

v _ (r) ji±vm + mm +m i nr) (51) 

The effective potential V mass i ve is plotted in Fig. 23 for varying values of the mass m 
of the scalar field. 
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Figure 23. The behavior of V mass i ve with the mass m of the scalar field. Here 
M = 1, q = 0.5, and 1 = 2. When the mass decreases, the height of the effective 
potential also decreases. 

The QNM's for the massive scalar field decay is computed using the WKB ap- 
proximation discussed in Section 4. They are given in the following figures, Fig. 24 
and Fig. 25. 
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Figure 24. The behavior of Re u with the mass of the scalar field m for M = 1, 
g = 0.76, and I = 1. Plots for n = and n = 1 are given. 
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Figure 25. The behavior of Im u with the mass of the scalar field m for M — 1, 
q = 0.76, and I = 1. Plots for n = and n — 1 are given. 

From Fig.24, the real part of the lowest QNM frequency increases with mass of 
the field. However, for the first overtone ( n — 1), the real part of QNM reaches a 
maximum and decreases rapidly with the mass of the field. Ohashi and Sakagami did 
an analysis on the QNM's for the Reissner-Nordstrom black hole with the massive 
scalar field in [49]. The real part of the u shows similar behavior as for the ur for 
n = for the Bardeen black hole. The imaginary part of the u decreases with the 
mass for the Bardeen black hole for both values of n. This is similar to the behavior 
shown for uij of the Reissner-Nordstrom black hole [49]. 

7 Summary 

We have studied the scalar perturbations of the Bardeen black hole. Quasinormal 
mode spectrum of the massless scalar field is computed for various values of the 
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charge q, mass M, and the spherical index I. The QNM spectrum is also computed 
for the Reissner- Nordstrom black hole by varying charge along with the Bardeen black 
hole and compared the behavior. 

We have also applied the unstable null geodesies of the black hole to compute the 
QNM frequencies in the eikonal limit ( I >> 1). Once again, a comparison is done 
with the QNM frequencies of the Reissner-Nordstrom black hole in the eikonal limit. 
We have also introduced the expansion method to compute u in inverse powers of L 
following the approach by [41]. 

Finally, the QNM frequency spectrum is computed for the massive scalar field. 
A discussion is also presented comparing the spectrum to the massive modes of the 
Reissner-Nordstrom black hole. 



8 Directions for further study 

There are several avenues to proceed from here to extend this work. For the Reissner- 
Nordstrom black hole, in the eikonal limit (/ ^> 1), the effective potentials for scalar, 
Dirac and gravitational perturbations are approximately the same. This is clearly 
presented in the paper by Ferrari and Mashhoon [39]. It would be interesting to do 
the analysis on gravitational perturbations and Dirac perturbations of the Bardeen 
black hole to see if similar behavior persists. 

Highly damped asymptotic QNM's of the black hole is not studied in this paper. 
As mentioned in the introduction, that is a very active area of research on QNM's. 
It would be interesting to extend this work to study u for large n values. Motl 
and Neitzke [50], did an analytical study to compute the asymptotic QNM's for the 
Reissner-Nordstrom black hole. It would be interesting to compare the asymptotic 
values of the Bardeen black hole with such values to understand how the nonlinear 
nature effects the physical properties of the black holes. 

It would be interesting to compute higher order expansion of uj in the eikonal limit 
using the approach by [41]. 

An extension of the Bardeen blackhole with a non-linear electromagnetic source 
presented by Ayon-Beato and Garcia [5] has been studied by Nomura and Tamaki 
[51]. The metric of this black hole is given by, 

ds 2 = -f(r)dt 2 + }\r)- l dr 2 + r 2 (d6 2 + sin 2 (9)d<p 2 ) (52) 

where, 

, n = , 2Mr 2 q 2 r 2 

Nomura and Tamaki concluded that the QNM frequencies become pure imaginary 
when n — > oo. It would be interesting to see if such behavior persists in all non-linear 
sources such as the one presented in this paper. 
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With regard to the QNM's for the massive scalar field, it would be interesting to 
do an analytical study if possible, to find a bound on the mass m when ojj becomes 
zero. 

Acknowledgements: SF like to thank R. A. Konoplya for providing the Math- 
ematica file for the WKB approximation. 
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